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Abstract 

We formulate a stochastic least-action principle for solutions of the incompressible Navier-Stokes 
equation, which formally reduces to Hamilton's principle for the incompressible Euler solutions in 
the case of zero viscosity. We use this principle to give a new derivation of a stochastic Kelvin 
Theorem for the Navier-Stokes equation, recently established by Constantin and Iyer, which shows 
that this stochastic conservation law arises from particle-relabelling symmetry of the action. We 
discuss issues of irreversibility, energy dissipation, and the inviscid limit of Navier-Stokes solutions 
in the framework of the stochastic variational principle. In particular, we discuss the connection of 
the stochastic Kelvin Theorem with our previous "martingale hypothesis" for fluid circulations in 
turbulent solutions of the incompressible Euler equations. 



1 Introduction 

Alternative formulations of standard equations can be very illuminating and can cast new light on 
old problems. As just one example, consider how Feynman's path-integral solution of the Schrodinger 
equation enabled intuitive new approaches to difficult problems with many-degrees-of-freedom, such 
as quantum electrodynamics and superfluid helium. In this same spirit, many different mathematical 
formulations have been developed for the equations of classical hydrodynamics, both ideal and non- 
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ideal. Recently, Constantin and Iyer [Tj have presented a very interesting representation of solutions of 
the incompressible Navier-Stokes equation by averaging over stochastic Lagrangian trajectories in the 
Weber formula [2 for incompressible Euler solutions. Their formulation is a nontrivial application of 
the method of stochastic characteristics, well-known in pure mathematics [3] (Chapter 6), in theoretical 
physics [HE] and in engineering modeling El [7]. The characterization of the Navier-Stokes solutions 
in pQ is through a nonlinear fixed-point problem, since the velocity field that results from the average 
over stochastic trajectories must be the same as that which advects the fluid particles. Constantin 
and Iyer have shown that their stochastic representation implies remarkable properties of Navier-Stokes 
solutions in close analogy to those of ideal Euler solutions, such as a stochastic Kelvin Theorem for fluid 
circulations and a stochastic Cauchy formula for the vorticity field. 

In this paper, we point out some further remarkable features of the stochastic Lagrangian formulation 
of [T] . Most importantly, we show that the nonlinear fixed-point problem that characterizes the Navier- 
Stokes solution is, in fact, a variational problem which generalizes the well-known Hamilton-Maupertuis 
least-action principle for incompressible Euler solutions [5]. We shall demonstrate this result by a 
formally exact calculation, at the level of rigor of theoretical physics. A more careful mathematical proof, 
with set-up of relevant function spaces, precise definitions of variational derivatives, etc. shall be given 
elsewhere. Closely related stochastic variational formulations of incompressible Navier-Stokes solutions 
have been developed recently by others [9] [TO] HJ and a detailed comparison with these approaches will 
also be made in future work. 

Our variational formulation sheds some new light on a basic proposition of pQ , the stochastic Kelvin 
Theorem which was established there for smooth Navier-Stokes solutions at any finite Reynolds number. 
We show that this result is a consequence of particle-relabelling symmetry of our stochastic action func- 
tional for Navier-Stokes solutions, in the same manner as the usual Kelvin Theorem arises from particle- 
relabelling symmetry of the standard action functional for Euler solutions [8]. This result strengthens 
the conjecture made by us in earlier work [121 113] that a "martingale property" of circulations should 
hold for generalized solutions of the incompressible Euler equations obtained in the zero- viscosity limit. 
Indeed, the stochastic variational principle for Navier-Stokes solutions considered in the present work 
is very closely similar to a stochastic least-action principle for generalized solutions of incompressible 
Euler equations that was developed by Brenier jTU [T5l [16] . One of the arguments advanced for the 
"martingale property" in [12] was particle-relabelling symmetry in a Brenier-type variational formula- 
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tion of generalized Euler solutions. That argument, however, did not distinguish an arrow of time, so 
that fluid circulations might satisfy the martingale property either forward or backward in time. It was 
subsequently argued in [13j that the backward-martingale property is the correct one, consistent with 
time-irreversibility in the limit of vanishing viscosity. The present work shows that a small but positive 
viscosity indeed selects the backward martingale property, as expected for a causal solution. 

2 The Action Principle 

The action principle formulated here for Navier-Stokes solutions involves stochastic flows 3J. The relevant 
flows are those which solve a backward Ito equation: 

i t x ro (a,i) = u ro (x ro (a,t),i)dt + V2vdW™(t),t < tt 

(1) 

x-(a,t / ) = a. 

V 

Here W ro (<), t 6 [ioji/] is a d-dimensional Brownian motion on a probability space (il,P,!F) which is 
adapted to a two-parameter filtration T\ of sub-cr-fields of J 7 , with to < t < t' < tf. Thus, W ro (s) — 
W ro (s') is T\ -measurable for all t < s < s' < t'. The constant v that appears in the amplitude of 
the white-noise term in the SDE fl} will turn out to be the kinematic viscosity in the Navier-Stokes 
equation. Note that, for such an additive noise as appears in (1), the (backward) Ito and Stratonovich 
equations are equivalent. 

In order to describe the space of flow maps which appear in the action principle, we must make a few 
slightly technical, preliminary remarks. The random velocity field u ro (r,£) in equation (fT]) is assumed 
to be smooth and, in particular, continuous in time, as well as adapted to the filtration T t , t < tf 
backward in time. It then follows from standard theorems (e.g. see Corollary 4.6.6 of [3]) that the solution 
x ro (a, t) of ([T]) is a backward semi-martingale of flows of diffeomorphisms. Conversely, any backward 
semi-martingale of flows of diffeomorphisms has a backward Stratonovich random infinitesimal generator 
F™(r, t) which is a spatially-smooth backward semi-martingale (e.g. see Theorem 4.4.4 of [3]). The class 
of such flows for which the martingale part of the generator is v^W" (t) and for which the bounded- 
variation part of the generator is absolutely-continuous with respect to dt coincides with the class of 
solutions of equations of form (1), for all possible choices of u ro (r, t) . Clearly, the random fields u CT (r, t) 
and x ro (a, t) uniquely determine each other. We consider here the incompressible case, where u ro (r, t) 
is divergence-free and x ro (a, t) is volume-preserving a.s. 
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The action is denned as a functional of the backward-adapted random velocity fields u ro (r,i) — or, 
equivalently, of the random flow maps x ro (a, t) — by the formula 

S[x] = J P{dvj)jjdt J d d r l -\u™{r,t)\ 2 (2) 

when this is well-defined and as +00 otherwise. The variational problem (VP) is to find a stationary 
point of this action such that x ro (a, tf) = a and x ro (a, to) = <p ra (a) for P — a.e. za, where y ro (a) is a 
given random field of volume-preserving diffeomorphisms of the flow domain. It is interesting that this 
problem is very similar to that considered by Brenier [14| I15[ 116] for generalized Euler solutions. The 
above problem leads instead to the incompressible Navier-Stokes equation, in the following precise sense: 

Proposition 1 . A stochastic flow x ro (a, t) which satisfies both the initial and final conditions is a 
solution of the above variational problem if and only i/u ro (r,£) solves the incompressible Navier-Stokes 
equation with viscosity v > 

d t u m + (u ro -V)u ro = -Vp m + v Au°, P-a.s. (3) 

where kinematic pressure p™ is chosen so that V-u ro = 0. 

Proof: Making a variation Su™ (r, t) in the random velocity field, the equation ([TJ) becomes 
cWx ro (a,£) = [5x ro (a,t)-V r u ro (x ro ,t) + <5u ro (x ro , i)]eft, t < tt 

(4) 

5x ro (a,£/) = 0. 

Since the VP requires that x ro (a, to) — (f m (a), one can only consider variations such that, also, 
<5x CT (a, t ) — 0. (We shall consider below an alternative approach with a Lagrange multiplier that 
permits unconstrained variations.) This equation may also be written as 

d t 6x™(a.,t) - (V r u ro (x ro ,i)) T <5x ro (a,<)di = 5u™(x™,t)dt (5) 

for t < tf and then easily solved by Duhamel's formula (backward in time) to give Skl^ (a, t) in terms 
of <5u ro (r,£). Since the martingale term vanished under variation, the process 5x ro (a, t) is of bounded 
variation and clearly adapted to the backward filtration T t s , t < tf. Conversely, any such flow variation 
will determine the corresponding velocity variation 5u ro (r, t) by the equation ([5]) directly. Lastly, note 
that the volume-preserving condition det(V Q x ra (a, t)) = 1 becomes 

V r -6x m (a™,t) = (6) 
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under variation, where a CT (r,£) is the "back-to-labels" map inverse to the flow map x ra (a, t). Be- 
cause these maps are diffeomorphisms, we see that the Eulerian variation of the flow map, 5x ro (r, t) = 
<5x tu (a CT (r, t), t), is an arbitrary divergence-free field. 

With these preparations, we obtain for the variation of the action j2]): 



SS[x] = [ P(dzu) [ 1 dt ( d d r u m {r,t)-6u™(r,t) 
J J to J 



P{dw) / d d a / u ro (x ro (a,i),t)- d t 5x ro (a, t) - <5x ro (a, £)-V r u ro (x ro , t)dt 

to 



P{dw) / d d a 





dt 







•<5x CT (a,t) 



(7) 



In the second line we employed (JSJ) . In the third line we integrated by parts, using the facts that 
<5x ra (a, tf) — (5x ro (a, to) = and that Jx ro (a, t) is a bounded variation process, so that the quadratic 
variation vanishes: dt(u ro (x ro , £), i5x CT (a, t)) = 0. We note that the final gradient term vanishes, because 
(5x ro (r,t) is divergence-free. We can evaluate the remaining term using the chain rule 



d t u ro (x ro ,i) = d t u m (x™,t)dt+ (x ro (a,o^)-V)u ro (x ro ,i), 



(8) 



in terms of the backward Stratonovich differential. This result can also be written using Ito calculus. 
Calculating from (T]| and © the quadratic variation 



/ 2vd t {WJ , {t), a Sj u ro (x w ,t)) = 2v A u ro (x ra , t)dt, 

one obtains the backward Ito equation 

d t u ro (x ro ,i) = [5 t u ro + (u ro -V)u ro -iyAu ro ](x ro ,t)di + \/2ri(dW ro (i)-V)u ro (x ro ,i), (9) 

A crucial point is that the martingale part of © vanishes when the expression is substituted back into 
(J7J), since both (Vu ro )(x ro ,i) and <5x ro (a, t) are adapted to the backward filtration T\ s , t < tf. Thus, 
the final result is 



£S[x] = / P(dn7) [ * dt [ d d r[9 t u ro + (u ro -V)u ro -i/Au ro ](r,t)^x ro (r,i). 



(10) 



Since the integrands are smooth in space and continuous in time and since the flow variation is an 
arbitrary divergence-free field, the theorem statement follows. □ 



There are alternative formulations of the VP which should be mentioned. Rather than performing the 
variation with the constraint x ro (a, to) = y ro (a), one can instead modify the action with a Lagrange 



multiplier term: 

S"[x,v ]=5[x]+ J P{dw) J d d a <(a).[x ro (a,i )-¥> ro (a)]. (11) 

Varying with respect to the Lagrange multiplier Vg 7 yields the constraint, whereas an unconstrained 
variation with respect to x CT yields 

5S'[x, v ] = J P{dw) J ' dtj d d r [d t u™ + (u CT - V)u CT - v A u ro ] (r, i)-c5x ro (r, i) 

+ y P{dm) J d d a K(a) -u ro (x ro (a,i ),t )]^x ro (a,t ). (12) 

The second term on the righthand side arises partly from the Lagrange multiplier term and partly from 
integration-by-parts in time. It follows that v J7 (a) can be identified as the Lagrangian fluid velocity at 
the initial time to and, likewise, ug^r) = v n (((^ ro ) _1 (r)) is the Eulerian fluid velocity at time to. Another 
alternative is to add a Lagrange multiplier term / P{dw) y t s dt J d d a 7r ro (a,i) lndet(V Q x ra (a,i)) for the 
incompressibility constraint and to allow variations over flows which are not volume-preserving. In that 
case 7r CT (a, t) is the Lagrangian pressure field and p ro (r, t) = 7r ra (a ra (r, t), t) is the Eulerian pressure. 

There are several mathematical questions that deserve to be pursued. Some technical issues remain, 
e.g. the precise degree of smoothness of solutions required to make the above argument fully rigorous, 
etc. It would also be very interesting to know under what conditions the solution of the VP corresponds 
to a minimum of the action and not just a stationary point. Although we have characterized the solutions 
of the VP, we have not proved cither their existence or their uniqueness. We just remark on the latter 
point that a unique stationary point certainly exists if the initial velocity Uo (r) is deterministic and if 
the Navier-Stokes equation has a unique solution u(r, t) over the time interval \to,tf] for that initial 
datum. This will be the case, for example, if the initial velocity is smooth enough and the Reynolds 
number Re = UL/v is low enough. In that case, the solution of the VP is also deterministic and is given 
by the corresponding Navier-Stokes solution. 



3 The Stochastic Kelvin Theorem 

We now mention a closely related result of [Tj: 

Proposition 2 (Constantin & Iyer, 2008) . The following two properties for a divergence- free 
velocity field u(r, t) are equivalent: (i) For all closed, rectifiable loops C and for any pair of times 
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to < t < a < t f , 

f u(a,t')-da= / P(dm) f u(r,t)-dr , (13) 
JC J /x™ t (C) 

where xj? t (a) are i/ie stochastic backward flows which solve equation (QP wf/j velocity u(r, t) /or times 
t < t' with final condition x£7 t ,(a) = a; and, (ii) the velocity u(r, t) satisfies the incompressible Navier- 
Stokes equation over the time-interval [to,tf]. 

This is just a slight restatement of Theorem 2.2 and Proposition 2.9 of pQ. The result (fT3")) is a stochastic 
version of the Kelvin Theorem on conservation of circulations for incompressible Euler solutions. Al- 
though circulations are not conserved for Navier-Stokes solutions in the usual sense, (fl~3|) states that 
circulations on loops advected by the stochastic Lagrangian flow are a martingale backward in time. 
This property of the Navier-Stokes solutions is closely related to the "martingale conjecture" of Eyink 
|12j for generalized Euler solutions obtained in the limit v — > 0. This connection will be discussed in 
detail in the next section. 

It is well-known that the Kelvin Theorem for incompressible Euler equations can be derived by 
the least-action principle as a consequence of an infinite-dimensional symmetry [5], called "particle- 
relabelling symmetry" and corresponding to the group of volume-preserving diffeomorphisms of the 
flow domain. This may be done by applying the general Noether Theorem relating symmetries and 
conservation laws. In this section, we shall show that the result of Proposition 2 can be similarly derived 
from the stochastic action-principle of section 1 as a consequence of particle-relabelling symmetry. We 
shall not make use of the Noether Theorem but, following Salmon [8], shall instead employ a more direct 
method of Lanczos |17) based on time-dependent symmetry transformations. 

Suppose given a smooth 1-parameter family {<p(a, t), t G [io,£/]} of volume-preserving diffeomor- 
phisms satisfying ip(a.,tf) = tp(a.,to) — a. Then any incompressible flow x(a, t) may be deformed into 
another such flow 

x v (a,i) = x(<p(a,i),i) (14) 
with initial and final values the same. It follows furthermore from (|14|) by chain rule that 



d t x v (a, t) = dfx(a, t) + (£(a, t)-V a )x(a, t)dt, (15) 

for a = tp(a,t). If x ro (a, t) is the solution of the stochastic equation |T]), then (|15jl implies that x?(a, t) 
also satisfies (fTJ) for the modified velocity field 

u?(r,t) - u ro (r,i) + (^(a ro ,t).V a )x ro (a ro ,i), (16) 



where we employ the shorthands a ro = a ro (r,£) and a ro = <£>(a ro (r, t), t). It is easy to see from (fTB]) 
that u?(r, i) is adapted to the backward filtration ^ 7 t / , t <tf whenever the original velocity u ro (r, t) is 
adapted. 

In infinitesimal form, cp(a,t) = a + eg(a, t) + 0(e 2 ) with V a -g(a, t) = and g(a, tf) = g(a, to) = 0. 
The formula (|16jl then yields the velocity variation 



<Su ro (r, t) = e(g(a ro , i)-V a )x ro (a ro , t) + 0(e 2 ). (17) 

The corresponding variation of the action is thus [see the first line of ([7|]: 

= <SS[x]=e y P{dvj) JJ dt J Ag(a,().w t, (a,()+0(e 2 ), (18) 

where w ro (a, t) is the stochastic Weber velocity [H[2j[8] 

w ro (a,t) = V a x ro (a,<)-u ro (x ro (a,i),i). (19) 

We can conclude that the 0(e) variation in (|18|) must vanish for any smooth divergence-free function 
g(a, t) with g(a, tf) = g(a, to) = 0. Taking limits of such functions, one may approximate a divergence- 
free distribution of the form 

g c , i)t < (a, r) = xit'.t] (r) £ 5 d (a - a') da' (20) 

for any closed, rectifiable loop C and any tf > t' > t > to- Here X[t,t']( T ) is the characteristic function 
of the interval [t, t'] and S d (a — a') is the Dirac delta-distribution. If we use the property of the Weber 
velocity that 

i w n (a',t).(Ja'= i u ro (r,f)-dr, (21) 

it then follows that 

/ P{dw) <f u ro (r,i')-dr = [ P(dw) <f u ro (r,t)-dr, (22) 

for any tf > t' > t > to- Taking the limit t' — > tf allows us to identify the above constant average 
as f c u(r,tf)-dr with u(r, t) = J P(dzu)u™ (r, t), since x CT (a,i/) = a P — a.s. When the velocity field 
that solves the VP is deterministic, i.e. corresponds to a unique Navier-Stokes solution u(r,i), then this 
result gives the statement (fT5|) of Proposition 2 for the special case where t' = tf. However, the VP may 
be applied not only over the entire interval [t , tf], but over any subinterval [t, t'] as well and this yields 
the general case. 



4 Irreversibility and the Zero-Viscosity Limit 



At first sight, it is strange to obtain the dissipative Navier-Stokes equation from a principle of least- 
action, which ordinarily leads to time-reversible equations. There is no paradox, however, since an 
"arrow-of-time" is built into the stochastic action- principle of Section 1. We may say that this is a 
causal variational principle, since labels are assigned at the final time and variations are over prior 
histories. The VP may be recast instead to be anti-causal, with fluid particle labels assigned at the 
initial time to and with flow maps solving a forward Ito equation: 

d t x ro (a,i) = u zu (x cu (a,t),t)dt + \/2^<iW ro (i), t > t ^ 
x CT (a,t ) = a. 

The random velocity field u ro (r,i) must now be adapted to the forward filtration !FL, t > to. An exact 
analogue of Proposition 1 holds, but with the conclusion that the velocity must satisfy 

d t u™ + (u ro -V)u ro = -Vp m - v A u ro , P-a.s. (24) 

or the negative-viscosity Navier-Stokes equation. An analogue of Proposition 2 also holds, in the form 



f u(a,t')-d&= / P(dw) f u(r,i)-dr 

JC J J x" „ (C) 



t < t' < t < t f , (25) 

/x- t (C) 

with circulations at the present time given anti-causally as averages over future values. The process of 
circulations in this case is a forward martingale. The proofs of all of the above statements follow by 
straightforward modifications of the previous arguments for the causal case. 

Conservation of energy is another property of Hamiltonian systems derived from the least-action prin- 
ciple, as a consequence of time-translation invariance. For example, consider a standard incompressible 
Euler fluid, with the action functional 

S[*]= l -J\t Jd d r |u(r,i)| 2 (26) 

Following the procedure of Lanczos [17] , one considers an arbitrary increasing function r(t) on the 
interval [to,tf], with r(to) — io and r(tf) = tf, and defines a modified flow 

x T (a,i) =x(a,r(t)). (27) 

When x(a, t) = u(x(a, t),t), then x T (a, t) satisfies the analogous equation with 

u T (r,t)=f(t)u(r,r(t)). (28) 



In an infinitesimal form, r(t) = t + £<5(t) + 0(e 2 ) with <5(t/) = <5(to) = 0, corresponding to a time- 
translation by a time-dependent shift. The variation in the velocity resulting from l|28p is 5u(r, t) = 
(d/dt)[S(t)u(r,t)], which implies a variation of the action 

SS[x] = f ' dtS(t) f d d r \\u{Y.t)\ 2 . (29) 
J t J 2 

From the stationarity of the action it follows that kinetic energy E{t) = hj d d r |u(r, t)\ 2 is conserved. 

This argument is not valid, however, for the stochastic action principle of Section 1. Indeed, if x ro (a, t) 
solves the stochastic equation ([T]) for backward Lagrangian trajectories, then the time-reparameterized 
flow (f27|) satisfies 

d t x?(a,t) = u?(x?(a,t),t)dt + y/2vf{t) dW ra (f), (30) 

where (r, t) is given by the analogue of (|28|) and W ro (t) is a Brownian motion defined on the same 
probability space as W ro (i). This is a consequence of standard results on time-change in stochastic 
differential equations (e.g. see [T5], Ch.IV, section 7). We thus see that the reparameterization (f2"7| 
leads to a flow map which is outside the class obeying an equation of the form ([1]) and for which 
the stochastic action ([2]) is formally +oo. Thus the argument leading to energy-conservation based on 
time-translation invariance of the action is no longer valid. 

One interest of the characterization of Navier-Stokes solutions via an action principle is that it may 
give some hint as to the character of their zero- viscosity limit. It was long ago conjectured by Onsager [19] 
that singular solutions of the Euler equations may result from that limit, relevant to the description of 
turbulent energy dissipation at high Reynolds numbers. For recent reviews, see [20l[2T]. The variational 
principle formulated for the Navier-Stokes solutions in the present work is similar to that of Brenier 
[16| for generalized Euler solutions, in which deterministic Lagrangian trajectories are also replaced by 
distributions over histories. There is other evidence to suggest that this may be a physical feature of 
the zero- viscosity limits of Navier-Stokes solutions, based upon recent results in a simpler problem, the 
Kraichnan model of random advection by a rough velocity field that is white-noise in time [221 [5] . Unlike 
the smooth- velocity case considered in [3] , it has been shown for the case of rough velocities that the 
solutions of the stochastic equations fl} and |23|) for backward and forward Lagrangian trajectories do 
not become deterministic in the limit as v — > |23l EH [25l [26l [27l [28l [29]. Instead, there are unique 
and nontrivial probability distributions on Lagrangian histories in the limit, a property referred to as 
"spontaneous stochasticity." This property is a direct consequence of Richardson's law of 2-particle 
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turbulent diffusion [30] and thus extends very plausibly to Navier-Stokes turbulence in the limit of large 
Reynolds numbers. 

The results discussed above helped to motivate the conjecture in [121 113] that a martingale property 
of circulations should hold for Eulcr solutions obtained in the limit v — > 0. The derivation of the 
stochastic Kelvin theorem in the present paper based on particle-relabelling symmetry is closely related 
to a similar argument in |12j for the martingale property of circulations in generalized Euler solutions 
(see section 4 there). Note, however, that it was incorrectly proposed in [T2] that circulations for such 
Euler solutions should be martingales forward in time, and this conjecture was only later emended to a 
backward-martingale property in [13] . The present work shows that the backward-martingale property 
is indeed the natural one, which could be expected to hold for dissipative Euler solutions obtained as 
the zero- viscosity limit of Navier-Solutions solutions. 

5 Final Remarks 

The stochastic Lagrangian representation of Constantin and Iyer pQ and our closely related variational 
formulation should clearly extend to a wide class of Hamiltonian fluid-mechanical models with added 
Laplacian dissipation. In a forthcoming paper [31] we prove the analogous results for several non- 
ideal (resistive and viscous) plasma models, including the two-fluid model of electron-ion plasmas, Hall 
magnetohydrodynamics (MHD), and standard MHD. As we shall show there, those models possess two 
stochastic Lagrangian conservation laws, one corresponding to the Alfven Theorem on conservation of 
magnetic flux [32] and another corresponding to a generalized Kelvin Theorem [331 EH] ■ hi following 
work we shall apply these results to important physical problems of magnetic reconnection and magnetic 
dynamo, especially in turbulent MHD regimes. 

As noted earlier, similar stochastic least-action principles have recently been proposed for the in- 
compressible Navier-Stokes equations [HI EH ELT] • A detailed discussion of the relation of these different 
variational principles to ours, as well as a rigorous treatment of the latter, will be the subject of future 
work. It is worth remarking that there is another variational principle for fluid equations, Onsager's 
principle of least dissipation 35, 36J, which determines the probability of molecular fluctuations away 
from hydrodynamic behavior in terms of the dissipation required to produce them. A modern formu- 
lation is presented in [37] . and [38) gives a rigorous derivation of Onsager's principle for incompressible 
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Navier-Stokes in a microscopic lattice-gas model. It would be interesting to know if any relation exists 
between the least-action and least-dissipation principles. 
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